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ABSTRACT 



This report presents an algorithm and HP-41CV programs for 
position determination with Loran-C chains. Additional 
computational routines include the ability to calibrate Loran 
station triplet data to a known benchmark and ITD's (Indicated 
Time Delay's), predict ITD's at given positions, compute the 
geodesic (similiar to great circle) bearing and distance from a 
fix to the destination and to compute the geodesic bearing and 
distance from any one location to another. Utility routines allow 
the user to transfer station pair data between the HP-41CV and 
magnetic cards, magnetic tape and an extended function/ memory 
module. 



A. Introduction 

The Loran system is a radio aid to navigation which utilizes 
the principle of hyperbolic fixing. The locus of points for which 
the difference in arrival time of synchronized signals from a pair 
of transmitters is constant determines a hyperbolic line of 
positions. The intersection of two hyperbolic lines of position 
from two pairs of stations determines position or a hyperbolic 
fix. That two pairs of stations are required for a fix does not 
necessarily mean that there are four separate stations, for one 
station of one pair may be colocated with one station of the other 
pair forming a Loran triplet (Figure 1). Triplets may be joined 
"end-to-end" by station colocation to form a Loran chain (Figure 
2). Loran chains are common on both the East and West coasts of 
the North American continent. 

The present day Loran-C operates at 100-kHz and is in use in 
the Atlantic, Pacific and Mediterranean areas. The computational 
algorithm and programs described herein can be used for position 
determination with Loran-C triplets. Further information on the 
history, development and operation of the Loran systems may be 
found in References 1 and 2. 



1 



M 



1 ' 



M 



2 



(a) Colocated Master Stations 




(b) Colocated Slave Stations 




(c) Colocated Master and Slave 



Figure 1. Loran Triplets. 
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Figure 2. 



Loran Chain of Five Loran Triplets. 
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B. User Routines 



1. Computational Routines; 

FI - The EXxing routine is the main program for calculating 
a Loran-C fix from indicated time delays. 

AS - The Alternate Solution routine will allow the second 

Loran fix solution to be computed. This routine toggles 
Flag 3 so that on subsequent fixes the FI routine will 
calculate the alternate solution. 

DN - The nestiMation routine stores the latitude and 
longitude of a fixed destination. 

HD - Computes the Heading and distance from the current fix 
to the destination stored by the DN routine. 

2. Manual Mode Routines: 

MI - Manual Input allows station data to be input and stored 
via the calculator keyboard. 

ED - Hcho Data is a utility routine for validating station 
triplet information stored in the calculator. 

3. Card Reader Routines; 

£S - Hard Store records station data onto magnetic cards. 

CR - Hard Head inputs station data from magnetic cards. 

CE - Hard £cho is a utility routine for validating station 
information stored on data cards. 

4. Extended Memory Routines; 

MS - XMEM Store records station data onto the extended 
memory module. 

XH - XMEM Hecall inputs station data from the extended 
memory module. 
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XU - XMEM Celete erases station data from the extended 
memory module. 

5. Tape Cassette Routines; 

IS. - Tape Store records station data onto the tape cassette. 

TE - Tape Eecall inputs station data from the tape cassette. 

TD - Tape Delete erases station data from the tape cassette. 

6. Utility Routines; 

DH - This routine is similiar to the HD option except that 

it computes the heading and distance from any origin to 
any destination. 

££ - PR is used to £Eedict the station ITD's that will be 
received at a given latitude and longitude. 

CA - CA is the Calibration option. Given the latitude and 
longitude of a known position and the indicated time 
delays from a Loran-C triplet, the stored station data 
are modified so that the FI routine (or AS) will 
compute the known position from the same time delays. 

SW - Dliitch data swaps the data of the two Loran stations 
stored in memory. 

Note; There are no specific routines that relate to HP-41CV 
printer operations. However, all input and output will be 
recorded on a printer if one is attached. 
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C. Recommended Kay Assignments 



It is recommended that the following HP-41CV user key 
assignments be made and recorded onto the program cards if the 
program is first prepared manually; 



CA ED 











CS 




CR 




CE 






TD 


MI 


SW 










TS 




TR 




DN 






XD 


PR 


DH 










XS 




XR 




HD 


AS 








FI 















The £I function is placed on the ENTER key as a reminder 
ENTER should not be used for data entry. To make the user key 
assignments on the HP-41CV, refer to the ASN function in the 
HP-41CV Owner’s Handbook and Programming Guide . 



that 



6 



D. User Instructions with Examples 
1. Manual I nput Ro utine - Ml 

This routine can be used to enter station data (Appendix 
C) to prepare data to be transferred to cards, extended memory, 
or tape cassette. It can also be used to enter station data 
manually if alternate data storage media are not available. 

Example using 9940W and 9940Y. (Note: The notation <9>, 
for example, means that you must press the gold shift key and 
then the 9 key.) 



See 


Key in 


Press 


1. 


XEQ"MI" 




2. ID? 


<9><9><4><0>W 


R/S 


3. CODE DELAY? 


11000 


R/S 


4. MS LAT? 


39.330662 


R/S 


5. MS LON? 


118.495637 


R/S 


6. SS LAT? 


47.034799 


R/S 


7. SS LON? 


119.443953 


R/S 


If desired, these 


station parameters can now be 


transferred to card, extended memory or 


tape using the CS, XS, 


or TS routines, respectively. Otherwise 


, repeat the steps above 


with the 9940Y station 
the remaining examples. 


data. These two 


stations will be used in 


S££. 


Key in 




1. 




XEQ"MI" 


2. ID? 


<9><9><4><0>Y 


R/S 


3. CODE DELAY? 


40000 


R/S 


4. MS LAT? 


39.330662 


R/S 
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5. MS LON? 



118.495637 



R/S 

6. SS LAT? 35.191818 R/S 

7. SS LON? 114.481743 R/S 

8. NEXT OPTION? 

At this point the station data for 9940W and 9940Y are 
stored in the calculator. The 9940Y data can be transferred to 
card, extended memory or tape cassette using the CS, XS or TS 
routines, respectively. 

[Advanced User Note: The first action taken in the MI 
routine is to XEQ"SW". SW is the data swap routine which 
exchanges to content of R16 - R25 with R26 - R35 . The incoming 
data are then stored in the R16 - R25 registers. The CS, XS and 
TS routines transfer the content of R16 - R25 to card, extended 
memory or tape cassette, respectively. If desired, the content 
of these registers can be swapped once more by using the SW 
utility routine.] 
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2. Echo Data Routine - ED 



This routine allows the user to review the station data 
resident in the calculator. 

Example: Load the station pairs 9940W and 9940Y using 
either the MI, CR, XR or TR routines. 



Kgy in Press 

1. XEQ"ED" 

2. ID: 9940W R/S 

3. CD: 11000 R/S 

4. MLT: 39.330662 R/S 

5. MLN: 118.495637 R/S 

6. SLT: 47.034799 R/S 

7. SLN: 119.443953 R/S 

8. BL: 2796.903 R/S 

9. ID: 9940Y R/S 

10. CD: 40000 R/S 

11. MLT: 39.330662 R/S 

12. MLN: 118.495637 R/S 

13. SLT: 35.191818 R/S 

14. SLN: 114.481743 R/S 

15. BL: 1967.302 R/S 

16. NEXT OPTION? 



Notation: CD = coding delay, M = master, S = slave, LT = 
latitude, LN = longitude and BL is the station pair baseline 
plus the secondary phase correction in microseconds. 
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3a. Card Store Routine„r_CS 



With the card reader attached, station data (in R16 - 
R25), which has been input using the MI, XR or TR routine, can 
be transferred to magnetic card using the XEQ"CS" command or by 
pressing the appropriate user defined key. 

Example using 9940W. 



See 


Key in 


Press 




1. 


XEQ"PD" 






2. WRITE: 9940W 


(Pass a blank card 
reader. Label the 


through the 
card track 


card 

"9940W") 


3. NEXT OPTION? 









To proceed with the remaining examples it is recommended 
that you also prepare a card for the 9940Y station pair. 

3b. Card Read Routine - CR 

With the card reader attached, XEQ"CR" or press the 
appropriate user defined key. This routine can be used to input 
the data for two station pairs, which must form a triplet. 



Example using 


9940W and 9940Y. 










See 


Key in 






1. 




XEQ"CR" 








2. 


1ST CARD 


(Pass the data card 
the card reader.) 


for 


9940W 


through 


3. 


2ND CARD 


(Pass the data card 
the card reader.) 


for 


99440Y 


through 


4. 


NEXT OPTION? 









10 



3c. Card Echo Routine - CE 

With the card reader attached, XEQ"CE" or press the 
appropriate user defined key. This routine is used to validate 
the content of data cards against the table in Appendix C. 
Example using 9940W. 

See Key, in 

1. XEQ"CE" 



2. STA. CARD {Pass one side of a data card [9940W] 

through the card reader) . 



ID; 9940W 


R/S 


CD; 11000 


R/S 


MLT; 39.330662 


R/S 


MLN; 118.495637 


R/S 


SLT; 47.034799 


R/S 


SLN: 119.443953 


R/S 


BL: 2796.903 


R/S 


NEXT OPTION? 





Notation: CD = coding delay, M = master, S = slave, LT = 
latitude, LN = longitude and BL is the station pair baseline 
plus the secondary phase correction in microseconds. 
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4a. Store Data in Extended_MemQrv =..XS 



With the extended memory module in the HP-41CV, station 
data (in R16 - R25), which has been input using the MI, CR, XR 
or TR routine, can be transferred to the module using the 
XEQ"XS" command or by pressing the appropriate user defined 
key. Example using 9940W. 

See Key in Press 

1. XEQ"XS" 

2. NEXT OPTION? 

Should the station pair already be in extended memory, the 
message DUP FL (duplicate file) will be displayed. If needed, 
the duplicate file may be erased using the XD routine. 

To proceed with the remaining examples it is recommended 
that you also store the 9940Y station pair. 

4b. Recall Data from Extended Memory - XR 

With the extended memory module installed, XEQ"XR" or 
press the appropriate user defined key to input the data for a 
station pair. 

Example using 9940W. (Note: The notation <9>, for example, 
means that you must press the gold shift key and then the 9 
key. ) 





See 


Key in 


Press 


1. 




XEQ"XR" 




2. 


ID? 


<9><9><4><0>W 


R/S 


3. 


NEXT OPTION? 
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4c. Delete Data from Extended Memory - XD 

XEQ"XD" or press the appropriate user defined key to 
delete the specific station pair data from the extended memory 
module. 

Example using 9940W. 

Kay_J.li Press 

1. XEQ"XD" 

2. ID? <9><9><4><0>W R/S 

3. NEXT OPTION? 

The message FL NOT FOUND will be displayed if the file you 
wish to delete in not in the extended memory. 

Note : The extended functions/ memory module will 
accommodate the data for 11 station pairs. The extended memory 
module will accommodate the data for an additional 22 station 
pairs. 
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5a. Store Data in the Tape Cassette. =._IS 



With the tape cassette attached to the HP-41CV, station 
data (in Rl6 - R25), which has been input using the MI, CR, XR 
or TR routine, can be transferred to the tape using the XEQ"TS" 
command or by pressing the appropriate user defined key. 

Example using 9940W. 

S££. Key in Press 

1. XEQ"TS" 

2. NEXT OPTION? 

Should the station pair already be in extended memory, the 
message DUP FL NAME (duplicate file name) will be displayed. If 
needed, the duplicate file may be erased using the TD routine. 

To proceed with the remaining examples it is recommended 
that you also store the 9940Y station pair. 

5b. Recall Data from the Tape Cassette -TR 

With the tape cassette attached, XEQ"TR" or press the 
appropriate user defined key to input the data for a station 
pair . 

Example using 9940W. (Note: The notation <9>, for example, 
means that you must press the gold shift key and then the 9 
key. ) 





See 


Key in 




1. 




XEQ"TR" 




2. 


ID? 


<9><9><4><0>W 


R/S 


3. 


NEXT OPTION? 
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5c. Delete Data from the Tape Cassette - TD 

XEQ"TD" or press the appropriate user defined key to 
delete the specific station pair data from the tape cassette. 
Example using 9940W. 

jS££ Ke y in Press 

1. XEQ"TD" 

2. ID? <9><9><4><0>W R/S 

3. NEXT OPTION? 

The message FL NOT FOUND will be displayed if the file to be 
deleted is not on the tape. 
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6. Loran-C Fixing Routines FI and AS 

Given the indicated time delay (ITD) from two station 
pairs which form a triplet, a Loran-C fix is obtained. 

Example: Load 9940W and 9940Y into the calculator using 
the MI, CR, XR or TR routine. The ITD on 9940W is 16019 
microseconds and the ITD of 9940Y is 42585 microseconds. Where 



are you? 





See 


K£y._in 


1. 






XEQ"FI 


2. 


ITD; 


9940W 


16019 


3. 


ITD: 


9940Y 


42585 


4. 


LAT; 


39.1419 




5. 


LON: 


115.5052 




6. 


NEXT 


OPTION? 





Since you are on a boat, you 
central Nevada at 39dl4'19" North 
Loran-C fix has two solutions, so 
alternate solution. 



2JL&&S. 

R/S 
R/S 
R/S 
R/S 

know that you cannot be in 
and 115d50'52" West. Every 
in this case you must use the 



Key in 

7. XEQ"AS" 

8. LAT: 35.0001 

9. LON: 125.0009 

10. NEXT OPTION? 



R/S 

R/S 



This is the proper solution at almost exactly 35 degrees 
North and 125 degrees West. Note that annunciator 3 (Flag 3) 
shows in the display indicating the alternate solution. If you 
should now repeat from Step 1, you will obtain the proper 
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solution iiranediately . 

The message "E: NO TRIPLET" will appear following Step 1 
if the data do not comprise a valid triplet. The latitudes and 
longitudes of each station pair at the vertex must agree 
exactly. Should this error occur, use the ED routine to review 
the resident station data. 

The message "E: ITD ERROR" will appear following Step 2 or 
3 indicating that the ITD you keyed in is inconsistent with the 
station parameters. Press R/S to be requeried for the ITD. 
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7a. Distance and Heading Routines DN and HD 



If you know the latitude and longitude of your 
destination, you may key these in and then see how far your fix 
is from your destination and what the geodesic heading 
(similiar to great circle heading) is to your destination. 

Example: Your destination is Moss Landing at about 36d48'N 
and 121d47'W. Your current fix is 35dN and 125dW (see the FI-SA 
example). First, key in your destination. 



See 


Key in 




1. 


XEQ"DN" 




2. DEST LAT? 


36.48 


R/S 


3. DEST LON? 


121.47 


R/S 


4. NEXT OPTION? 




R/S 



The destination is now stored in the calculator and will remain 
unchanged until you use either the DN or DH options. Also, the 
latest fix is stored and will remain unchanged until you use 
either the FI, AS or DH options. Now, find the distance and 
bearing from the latest fix (see the FI-AS example) to Moss 
Landing. 



See 



Key in 



Press 



1 . 



XEQ"HD" 



2. N.MI: 190.38 R/S 

3. BRG: 54.3411 R/S 

4. NEXT OPTION? R/S 

So the distance to Moss Landing is 190.38 nautical miles 
at a heading of 54d34'll". 
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7b. Distance and Heading Routine - DH 



Given the latitude and longitude of an origin and 
destination, this routine will find the distance and heading 
from one to the other. 

Example: How far, and in what direction, is Corvallis, 
Oregon (44d34'N, 123dl6 *W) from Cupertino, California (37dl9N, 
122d02 'W) ? 



See 


K£^.Id 






XEQ"DH 


ORIG 


LAT? 


37.19 


ORIG 


LON? 


122.02 


DEST 


LAT? 


44.34 


DEST 


LON? 


123.16 


N.MI 


: 438.32 




BRG: 


353.0259 




NEXT 


OPTION? 





R/S 

R/S 

R/S 

R/S 

R/S 

R/S 



Thus the distance 



is 438.32 nautical miles and the 



direction of Corvallis 



from Cupertino 



is 353d02'59". 
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8. ITD Prediction R outine - PR 

As an aid to identification, this routine will allow the 
user to determine what ITD's should be received at a given 
location. 

Example: Suppose that you know you are somewhere near 
latitude 35 North and longitude 125 West but are not sure what 
ITD's you should be receiving from 9940W and 9940Y. To 
determine these ITD's, proceed as follows: 



See 




Key in 


Press 






XEQ"PR'' 




LAT? 




35 


R/S 


LON? 




125 


R/S 


9940W: 


16019.35 




R/S 


9940Y: 


42584.71 




R/S 


NEXT OPTION? 







You should expect to receive an ITD of 16019.35 from 9940W 
and an ITD of 42584.71 from 9940Y. 
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9. Calibration Routine - CA 

This routine will allow the user to calibrate the Loran 
data in the calculator to a known position when the indicated 
time delay (ITD) is known for each station pair. 

Example: Suppose you are receiving an ITD of 16308 from 
9940W and 42800 from 9940Y. These ITD's would tell you that 
your location is 36d47'55"N and 121d47'll"W. However, you know 
that your position is bench marked to be at 36d47'36"N and 
121d46 ' 58 ”W, and you wish to calibrate your calculator so that 
the ITD's of 16308 and 42800 will give you the latter fix 
instead of the former. Proceed as follows: 





See 


K£^-.ln 


Press 


1. 




XEQ"CA" 




2. 


LAT? 


36.4736 


R/S 


3. 


LON? 


121.4658 


R/S 


4. 


ITD 9940W: 


16308 


R/S 


5. 


ITD 9940Y: 


42800 


R/S 


6. 


NEXT OPTION? 







Entering 16308 and 42800 into the FI routine will now give 
you a fix at 36d47'35"N and 121d46'59"W. The small discrepancy 
between this fix and the bench mark is due to assumptions made 
in the fixing algorithm. 

Calibration is achieved by modifying the Master/Slave 
baseline (BL in the CE and ED routines). See Section F. 
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10. Switch Data Registers .RoutiiiB-^-SH 



The SW utility allows the user to swap the station data 
stored in R16 - R25 with the data stored in R26 - R35 . Whenever 
the MI, XR or TR routine is used, the SW routine is invoked 
prior to the loading of the data; the incoming data are then 
placed in R16 - R25 . The consequence is that the first station 
pair data reside in R26 - R35 and the second station pair data 
reside in Rl6 - R25 . 

One user application of SW would be to change the order of 
the station ID query in the FI routine (this also affects the 
order of determination of the solution and alternate solution) . 
Another user application would be to output both resident 
station pairs to card, extended memory or tape cassette using 
the CS, XS or TS routines, respectively. To accomplish this, 
first use the CS, XS, or TS routine; then XEQ"SW" (note that 
the ID of the station data in R16 - R25 appears in the display 
instead of the NEXT OPTION? prompt) ; and finally use the CS, XS 
or TS routine once more. 
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11a. fiaCQ-Cdi nq . the Loran-C Program onto Magnetic Cards 



(1) Attach the card reader to the HP-41CV. 

(2) Place the calculator in the USER mode. 

(3) Press the PRGM key to place the calculator in the 
program mode. 

(4) Pass one side of a blank magnetic card through the 
card reader. Then follow the display prompts until 
nine program cards (17 tracks) have been recorded. 

(5) Press the PRGM key once more to leave the program 
mode. 

(6) To record a status card, XEQ"WSTS". Then pass a blank 
card through the card reader following the display 
prompts. 

11b. Recording the. Loran-C Program onto Magnetic Tape 

(1) Attach the tape cassette to the HP-41CV. 

(2) Place the calculator in the USER mode. 

(3) Press the alpha key, key in the word LORANC, press the 

alpha key once more, and then XEQ"WRTP". 

(4) To record the program status, press the alpha key, key 

in the word STATUS, press the alpha key once more, and 

then XEQ"WRTS". 
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12a. Loading the Loran-C Program from Magnetic Cards 



(1) Attach the card reader to the HP-41CV. 

(2) Clear program memory: Turn the calculator off, then, 
while pressing the left arrow (erase) key down, turn 
the calculator on. 

(3) Place the calculator in the USER mode. 

(4) Read in the STATUS card. The status card will set the 
calculator to SIZE 42. 

(5) Read in the nine program cards (17 tracks). 

12b. Loading the Loran-C program from Tape Cassette 

(1) Attach the tape cassette to the HP-41CV. 

(2) Clear program memory. (See Step 2 above). 

(3) Place the calculator in the USER mode. 

(4) Press the alpha key, key in the word STATUS, press the 
alpha key once more, and then XEQ"READS". The status 
file will set the calculator to SIZE 42. 

(5) Press the alpha key, key in the word LORANC, press the 
alpha key once more, and then XEQ"READP". 
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E. The Loran-C Fixing Algorithm 



The principles of Loran lines of position (LOP's) and 
fixing are adequately covered in Reference 1 and will not be 
repeated here. 

The basic Loran-C equation [Ref. 4] can be written as 
ITD = [Tg + p(Tg)] - [T^ + P(T^)] + [Tg + p (T^) ] + 6 (1) 

where 

ITD is the "indicated time difference" in microseconds, 

T^ is the distance, in microseconds, from the master 
to the receiver, 

Tg is the distance, in microseconds, from the slave to 
the receiver, 

T„ is the distance, in microseconds, between the master 

D 

and the slave, 

(5 is the assigned station pair coding delay, in micro- 
seconds, and 

p(T) is the secondary phase correction, in microseconds, 
for a surface seawater path of length T . 

The quantity 



At = Tg + p(Tg) + 6 

is a constant for each master/slave pair. The quantity Tg is 
computed from the positions of the master and slave using the 
reverse solution algorithm (Section H) at the time of manual 
data input (Routine MI) . 
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The following World Geodetic System 1972 (WGS 72) values 
have been adopted for Loran-C navigation [Ref. 4]; 

Vq = 299792458 raeters/second is the velocity of light 
in free space, 

n = 1.000338 is the index of refraction of the 

surface of the earth for standard atmosphere 
and 100 kHz electromagnetic waves, 
a^ = 6378135.000 meters is the equatorial radius of 
the earth 

and f = 1/298.26 is the flattening factor (1-b/a^, where 

b is the polar radius) of the earth. 

Accurate formulas for computing the secondary phase cor- 
rection p(T) are contained in Reference 4, but for use in the 



HP-41CV, 


the 


form 












p(T) = a^/T + a^ + 


a2T 


is used. 


where T 


is in microseconds 


and 


1. 


For 


1 V 


537 ysec: 








^0 = 


129, 








^1 = 


-0.408, and 








^2 = 


0.0006458. 




2. 


For 


T < 


537 ysec: 





a^^ = -0.011, and 
= 0.00033. 



26 



On the surface of a sphere, a hyperbolic line of position 



(LOP) can be represented by the equation [Ref. 1, page 175] 



tan r 



cos 2a - cos 2c 

sin 2c cos co + <; sin 2a 



( 3 ) 



where the origin of the coordinate system is at the prime focus 
of the spherical hyperbola, 2c is the spherical arc joining the 
foci, 2a is a constant for any one hyperbola, and r and w 
are the spherical coordinates of a point on the hyperbola. If 
the base line of the coordinate system is the arc joining the 
foci the to is the spherical polar angle from the baseline to 
a point P on the spherical hyperbola and r is spherical 
polar distance (or arc) from the prime focus to P . Using the 
Loran system we take ^ = +1 if the prime focus is at a master 
station and <; = -1 if the prime focus is at a slave station. 

If we let V = Vq/h be the velocity of lOOkHz electro- 
magnetic radiation at the earth's surface then, for a spherical 
earth, we can relate the parameters in Equations 1 and 3 as 
follows : 



2c = vTg/a^ , 

and 

2a ' v(Tg - T^)/ag . 

Using the spherical approximation for now, we see that 2c is 
known for any Loran pair. The "indicated time delay" ITD is 
measured by the receiver at point P , and to determine a 
hyperbolic line of position we must determine 2a, but 
cannot be computed from Equations 1 and 2. If a^ were zero 
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in Equation 2, then it would be possible to determine Tg - 
uniquely. As an approximation we use the following parameters 
in Equation 2: 



a^ = -0.321 , 

and ^2 — 0.000635 . 

These values have been obtained by setting a^ = 0 and determin- 
ing a^^ and a 2 by linear regression of the T > 537 values over 
the interval of 1000 < T < 8000. This approximation is quite 
good (within 0.03 ys) for distances up to 10,000 microseconds 
where small changes in the LOP ' s can cause large position errors. 
At short distances the error increases from 0.05 ys at 1000 ys 
to 0.58 ys at 10 ys; although these errors are large for small 
distances, the LOP ' s are not as sensitive to these changes as 
they would be at large distances. When this approximation is 
substituted into Equation 1, we obtain 

■ [’'m + ' 

or 

Tg - T^ = (ITD - At)/(1 + a2) (4) 

and hence 2a = v(T„-T.J/a is determined for use in the 

S M e 

spherical approximation. 

Consider a Loran-C triplet with the master stations co- 
located. Let and ^2 denote the azimuth angles of slave 

1 (Sj^) and slave 2 (S 2 ) / respectively, measured from North 
toward the East from the master stations (M) (see Figure 3) . 
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Further, let a and r be the azimuth and spherical polar 
arc (distance) of the receiver (R) from M . For this geometry 
Equation 3 can be written in the form 



where 



tan r . 
1 



B. 

1 

cos (a - ' 



A. = C- sin 2a. , 

11 1 ' 



( 5 ) 



and 



B. = cos 2a. - cos 2c. , 

1 1 1 ' 

C. = sin 2c. 

1 1 



for the i^^ Loran pair, i = 1,2. Since r^ = r 2 = r, we can 
eliminate tan r between the two equations. The resulting 
equation can be rewritten as 



where 



and 



COS a 


+ S sin 


P 

II 




= 


cos ^2 


- ^ 2^1 


^1 ' 


(N 

u 
1 — 1 
PQ 

II 


sin ^2 


- B 2 Cj^ sin 


^1 ' 


= 


- ^1^2 


• 





If we define p > 0 and y by the equations 



( 6 ) 



p cos y = C , (7) 

and p sin y = S , 

then 

p = /c^ + , 

and y = qatn (S,C) . 
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Figure 3. Geometry of a Loran Triplet and a Receiver. 
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Here the function qatn(y,x) is the arctangent of y/x adjusted 
for the proper quadrant according to the signs of x and y . 
A compact form of this function is 



qatn(y,x) = tan 



-1 



X + 10 ^ t (x = 0?) 



+ TTt (X < 0? ) 



where t(z) = 1 

and t(z) = 0 



when z is true 
when z is false, 



When convenient we will use the notation qatn(y/x) inter- 
changeably with qatn(y,x). Now we can substitute Eq. (7) into 
Eq. (6) and solve for 

a = Y ± cos ^(K/p) (8) 

to obtain the azimuth angle a of the two points of inter- 
section of the spherical hyperbolic LOP ' s . Finally we can 
obtain a value for r by substituting each a into Eq. (5) 
for either i = 1 or i = 2. We find that 



r 



qatn 



B. 

1 

C . cos (a-? . ) + A. 
1 1 1 



for i 



1 or 2 . 



(9) 



The distance and azimuth from M or the triplet vertex can be 
converted into the latitude and longitude of the two possible 
positions of r . 

The fixing algorithm then uses a and r in the direct 
solution algorithm of spheroidal geodesy (Section G) . 
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F. The Calibration and ITD Prediction Algorithms 

Calibration can be achieved when an ITD is measured at a 
known bench marked position. From the bench marked position 
and the known master and slave positions, the quantities 
Tj^ + P(T^) + P(Tg) can be computed using the reverse 

solution algorithm (Section H) and the accurate secondary phase 
correction formula (Eq. 2) . Equation 1 can then be solved for 
T^ + p(Tj,) to obtain a modified baseline. This modified base- 
line is stored and then used instead of the true baseline in 
subsequent computations. The affect on the accuracy of fixes 
using this modified baseline with positions far removed from 
the bench mark has not studied. 

The ITD prediction algorithm is a direct application of 
Equation 1. A known position, together with the known master 
and slave positions, is used to compute the quantities 
T^ + P(Tj^) and T^ + p(Tg). When these values, along with a 
computed or calibrated baseline, T„ + p(T„) , are substituted 

into Equation 1, a predicted ITD is obtained. 
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G. 



The Direct Solution Algorithm 



This direct solution algorithm is a modification of the 
second order in flattening (f) algorithm given by Thomas 
[Ref. 5, pp. 7-8]. Thomas' notation has been followed as 
closely as possible for ease of comparison of the algorithms. 

The qatn function is defined in a previous section. East 
longitudes and South latitudes are negative. We are given the 
point ^ the spheroid, where the geodetic 

latitude and longitude (geographic coordinates) ; the forward 
azimuth distance S to a second point P2^‘P2'^2^ ' 

from these we are to find the geographic coordinates (i> 2 > ^2 
the back azimuth 0 . 2 -^ • The given quantities are ^-^2 

and S . No assumptions about the relative location of 
and P 2 are required. The modified direct solution algorithm 
is : 

= tan ^[(l-f)tan 6j^J, M = - sin a ^^2 ®]_ > 

= fM, C 2 = f(l-M^)/4 , 

D = {I-C 2 ) (I-C 2 - Cj^M) , P = C2[l + (1/2)C^M]/D , 

N = cos cos a ^^2 r = qatn(N, sin 6^^) , 

d = S/(a^D) , u = 2{aj^-d) , W = 1-2P cos u , 

V = cos(u+d) , X = C 2 sin d cos d(2V^-l) , 

Y = 2PVW sin d,Aa=d+X-Y, 

n 2 1/2 

K = [M + (N cos-Ao sin 6^^ sin Aa) ] ' , 

©2 = tan (sin 6^^ cos Ao + N sin Aa)/K] , 

An = qatn(- sin Aa sin a ^^2 / ~ ^12^ ' 

H = Cj^(l-C 2 ) Aa - Cj^C 2 sin Aa cos(2aj^ - Aa) , 

A 2 = Aj + An - H , 
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= qatn[-M, - (N cos Aa - sin 6 ^^ sin Ao) ] , 

(J>2 = tan~^[tan 02/(l-f)] 

Details of the modifications made to Thomas' algorithm are 
contained in Reference 3 . The algorithm above has been further 
modified so that Eastern longitudes, rather than Western 
longitudes, are negative. 



34 



H . The Reverse Solution Algorithm 



This reverse solution algorithm is a modification of the 
second order in flattening (f) algorithm given by Thomas 
[Ref. 5, pp. 8-10]. Thomas' notation has been followed as 
closely as possible for ease of comparison of the algorithms. 
The qatn function is defined in a previous section. East 
longitudes (A) and South latitudes (<))) are negative. We are 
given the points (cj) , A , P2(<}»2/^2^ spheroid and are 

to find the distance S between the points and the forward 
and back azimuths, U 2 j^ . Given quantities are 

Aj^, (P 2 ^2 ' assumptions about the relative location 

of Pj^ and P 2 are required. The modified reverse solution 
algorithm is: 

0 ^ = tan~^[ ( 1 -f ) tan (J)^] , i = 1,2 , 

AA = A 2 - Aj^ , A0^ = (02 - 0 jl)/2 , + ^ 2 ^/^ > 

H = cos^ A0 - sin^ 0 , L = sin^ A0^ + H sin^ (AA/2) , 

mm m 

d = 2 sin ^/l , U = 2 sin^ 9 cos^ A 6 /(1-L) , 

m m 

V = 2 sin^ A0 cos^ 0/L,X=U+V,Y=U-V, 
m m 

T = d/sin d , D = 4T^ , E = 2 cos d , A = DE , 

C = T - (A-E)/2, n^ = X(A+CX) , 

B = 2D, n 2 = Y(B+Ey) , n^ = DXY , 

= f^(n^ - 03 ) /64 , 6 ^d = f(TX-Y)/4 , 

S/a^ = (T - 6 j^d + 62 ( 3 ) sin d , M = 32T - (20T-A)X - (B+4)Y , 

F = 2Y - E(4-X) , G = fT/2 + f^M/64 , Q = - (FG tan AA)/4 , 

AA ' = (AA + Q)/2 , 

m 
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qatn(sin A0 cos AX’ , cos 0 sin AA') , 

^ mm mm 

qatn(- cos A0 cos AA' , sin 0 sin AA') , 

^ mm mm 

^ ~ ^2 

Details of the modifications made to Thomas' algorithm 
are contained in Reference 3. The algorithm above has been 
further modified so that Eastern longitudes, rather than 
Western longitudes, are negative. 
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I. Program ...Ac.g u racy 



The direct and reverse solution algorithms are equivalent to 
the second order flattening algorithms given by Thomas (Ref. 5) ; 
the parameters of the WGS 1972 spheroid are used. The reverse 
solution algorithm reproduces the baselines provided by the 
Defense Mapping Agency for all 40 Loran-C stations to within 0.15 
meters (the average deviation is -0.031 meters, DMA minus HP-41CV, 
with a standard deviation of 0.037 meters) and to within 0.01 
microseconds, including the secondary phase correction for an all 
seawater path. The reverse solution algorithm is also used to 
generate predicted ITD's; these are presumed to be within 0.01 
microseconds, also. 

The fixing algorithm uses the direct solution algorithm with 
the azimuth and distance of the fix from the vertex of the Loran-C 
triplet computed from Equations 8 and 9 as inputs. Equations 8 and 
9 are based upon spherical geometry and include an approximation 
to the secondary phase correction for an all seawater path. The 
largest source of error is the assumption that the azimuth and 
distance to the fix are accurately represented by this spherical 
approximation. This approximation has not been rigorously tested, 
however it is possible to use the reverse solution algorithm to 
predict the ITD's that will be received at a given position and 
then to enter these ITD's into the fixing algorithm to determine 
how accurately the fixing algorithm reproduces the original 
position. The distance between the fix and the original position 
can be determined using the HD algorithm. Tables 1, 2 and 3 were 
producted in this manner. Similiar tables with different station 
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pairs are given in References 3, 6 and 7. It is felt that the 
results are accurate enough to not warrant the inclusion of an 
iterative improvement routine. Of the samples in Tables 1, 2 and 
3, the largest error, 0.22 n.mi., is the first entry in Table 1. 
From the chart LCNC-2, it estimated that the angle of intersection 
of the two hyperbolic lines of position is about 5 degrees and so 
an error of 0.22 n.mi. should not be unexpected. 
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Table 1. Colocated Master Stations 



Position 
Lat Long 


Predicted ITD's 
9940W 9940X 


HP- 4 lev Fix 
Lat(N) Long(W) 


Error 

n.mi 


31° 


123° 


16413.28 


27570.93 


30°59’47" 


123°00’03" 


0.22 


37° 


126° 


15610.11 


27020.50 


36°59'53" 


126°00'09" 


0.17 


42° 


129° 


13881.78 


27285.58 


42°00'00" 


129°00'00" 


0 


44° 


132° 


13180.89 


27371.19 


44°00'00" 


132°00'01" 


0.01 


00 

o 


135° 


12301.25 


27552.06 


48°00’01" 


135°00'03" 


0.03 


0 

O 


138° 


12068.67 


27584.22 


50°00’01" 


138°00'04" 


0.05 



Table 2. Colocated Slave Stations 



Position 
Lat Long 


Predicted ITD's 
9940W 5990Y 


HP-41CV Fix 
Lat(N) Long(W) 


Error 

n.mi 


31° 


123° 


16413.28 


27177.18 


31°00’03" 


122°59'58" 


0.06 


37° 


126° 


15610.11 


27403.20 


37°00'01" 


125°59'59" 


0.01 


0 

CM 


129° 


13881.78 


27955.45 


42°00'00" 


129°00'00" 


0 


44° 


132° 


13180.89 


28512.90 


44°00'00" 


132°00'00" 


0 


0 

00 


135° 


12301.25 


29413.61 


48°00'00" 


134°59'59" 


0.01 


50° 


138° 


12068.67 


29816.84 


50°00'00" 


137°59'59" 


0.02 
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Table 3. Colocated Master and Slave Stations 



Position 
Lat Long 


Predicted ITD’s 
5930Y 9960W 


HP-41CV Fix 
Lat(N) Long(W) 


Error 

n.mi 


44° 


o^ 

o 


29864.46 


11685.15 


44°00’00" 


63°00'00" 


0 


41° 


66° 


30585.61 


12946.91 


41°00'00" 


66°00’00" 


0 


39° 


69° 


31020.46 


14111.31 


39°00'00" 


69°00’00" 


0 


35° 


72° 


31064.57 


15139.48 


35°00’00" 


72°00'00" 


0 


30° 


75° 


31040.82 


15610.46 


29°59’59" 


75°00'01” 


0.02 


26° 


O 

CO 


31106.20 


15858.46 


25°59’57" 


78°00'02" 


0.05 
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Appendix 


A: Program Storage Allocations, Flag Usage 




and 


Program Listing 




Program Size: 42 






Registers ; 








ROO - R13: 


Scratch 


storage 




Rl4: Flattening 






2nd Sta. 


Pair 


Variable 


1st Sta. Pair 


R15 : 




2a 




R16 : 




ID 


:R26 


R17: 




Tb + p(Tb> 


:R2 7 


R18: 




2c 


:R28 


R19 : 




CD 


:R29 


R20 : 




0 master 


;R30 


R21: 




A master 


:R31 


R22 : 




a master-slave 


:R32 


R23: 




0 slave 


:R33 


R24 : 




A slave 


:R34 


R25: 




a slave-master 


:R35 






2a 


:R36 


R37: 




0 fix 




R38: 




A fix 




R39 : 




a fix-dest 




R4 0 : 




0 dest 




R41 : 




A dest 
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Flag Usage : 



FOO : 

FOl: 
F03: 
F05: 
F06 : 
F0 7: 
F08: 
F0 9 : 
FIO : 

F14: 



. Requery erroneous ITD at LBL06 
. Echo two data sets at LBL ED 
PR and CA interlock 
Set for 2nd solution 

Set if slave at vertex of 2nd station pair 
Set if slave at vertex of 1st station pair 
Set for XMEM functions, clear for tape functions 
DN and DH interlock 
Set if DMS conversion required 
• ^ interlock 

. Input interlock in FI vertex check 
PR and ^ loop control 
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Program Listing 



D 

< 

C 

•H 

C/5 



D 

< 
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a 

IS 



CD 

c 

•H 

C/3 
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IS 
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Append ix B; station coverage 



Station No. of Pai 



4990 2 

5930 2 

5990 3 

7930 3 

7960 2 

7970 4 

7980 4 

7990 3 

8970 3 

9940 3 

9960 4 

9970 4 

9990 3 



s. Location 

Central Pacific 
East Coast, Canada 
West Coast, Canada 
North Atlantic 
Gulf of Alaska 
Norwegian Sea 
Southeast U.S.A. 
Mediterranean Sea 
Great Lakes 
West coast, U.S.A. 
Northeast U.S.A. 
Northwest Pacific 
North Pacific 
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Appendix C; LORAN-C STATION DATA 



The following list contains the pertinent parameters for each 
Loran-C station pair. This list was compiled from the data in 
Reference 4. Each column contains the following information: 

1. The Loran-C station pair designator. 

2. The coding delay. 

3. The master station latitude. 

4. The master station longitude. 

5. The slave station latitude. 

6. The slave station longitude. 



In this list, positive longitudes are West, negative 
longitudes are East, positive latitudes are North and negative 
latitudes are South. In columns 3 through 6 the latitudes and 
longitudes appear to be in decimal form, but the actual format is 
DDD.MMSSFF (which is compatible with the HP-41CV D.MS or H.MS 
input mode) where 



DDD designates degrees, 

MM designates minutes, 

SS designates seconds, and 
FF designates hundredths of seconds. 
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Appendix D; CQLQCAI£D-£TATIQNS 



5930 MASTER 


— 


9960W 


5930X 


- 


9960X 


5930Y 


- 


7930Z 


5990X 


- 


7960Y 


5990Y 


- 


9940W 


7930W 
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7970Y 


7930X 
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7970 MASTER 


7960X 
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9990Z 


7980Z 
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9960Y 


7980 MASTER 


- 


8970W 


8970X 


- 


9960 MASTER 


8970 MASTER 


_ 


9960Z 



Slave stations are denoted with a letter 
suffix. Master stations are so designated. 
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